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Abstract 

Three-jet production in deep inelastic ep scattering and photoproduction was 
investigated with the ZEUS detector at HERA using an integrated luminosity of 
127 pb _1 . Measurements of differential cross sections are presented as functions 
of angular correlations between the three jets in the final state and the proton- 
beam direction. These correlations provide a stringent test of perturbative QCD 
and show sensitivity to the contributions from different colour configurations. 
Fixed-order perturbative QCD calculations assuming the values of the colour 
factors Cf-i Ca and Tp as derived from a variety of gauge groups were com- 
pared to the measurements to study the underlying gauge group symmetry. The 
measured angular correlations in the deep inelastic ep scattering and photopro- 
duction regimes are consistent with the admixture of colour configurations as 
predicted by SU(3) and disfavour other symmetry groups, such as SU(iV) in the 
limit of large N. 
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1 Introduction 



Quantum chromodynamics (QCD) is based on the non-Abelian group SU(3) which induces 
the self-coupling of the gluons. Investigations of the triple-gluon vertex (TGV) were 
carried out at LEP [1,2] using angular correlations in four-jet events from Z° hadronic 
decays. At HERA, the effects of the different colour configurations arising from the 
underlying gauge structure can be studied in a clean way in three-jet production in neutral 
current (NC) deep inelastic scattering (DIS) and photoproduction (jp). 

Neutral current DIS at high Q 2 (Q 2 ^> Aq CD , where Q 2 is the virtuality of the exchanged 
photon) up to leading order (LO) in the strong coupling constant, a s , proceeds as in 
the quark-parton model (Vq — > q, where V = 7* or Z°) or via the boson-gluon fusion 
(Vg — > qq) and QCD-Compton (Vq — > qg) processes. Photoproduction is studied at 
HERA by means of ep scattering at low four- momentum transfers (Q 2 ~ 0). In 7p 
reactions, two types of QCD processes contribute to jet production at LO [3,4]: either 
the photon interacts directly with a parton in the proton (the direct process) or the photon 
acts as a source of partons which scatter off those in the proton (the resolved process). 

A subset of resolved subprocesses with two jets in the final state are described by dia- 
grams with a TGV; however, such events are difficult to distinguish from two-jet events 
without such a contribution. Three-jet final states in direct jp processes also contain 
contributions from TGVs and are easier to identify. Since three-jet production in NC DIS 
proceeds via the same diagrams as in direct 7/?, such processes can also be used to inves- 
tigate the underlying gauge symmetry. Examples of diagrams contributing to four colour 
configurations are shown in Fig. 1: (A) double-gluon bremsstrahlung from a quark line, 
(B) the splitting of a virtual gluon into a pair of final-state gluons, (C) the production 
of a qq pair through the exchange of a virtual gluon emitted by an incoming quark, and 
(D) the production of a qq pair through the exchange of a virtual gluon arising from the 
splitting of an incoming gluon. 

Other possible diagrams and interferences correspond to one of the four configurations. 
The production rate of all contributions is proportional to the so-called colour factors, 
Cf, Ca and Tp, which are a physical manifestation of the underlying group structure. 
For QCD, these factors represent the relative strengths of the processes q — > qg, g — > gg 
and g — > qq. The contributions of the diagrams of Fig. 1 are proportional to C F , CfCa, 
CpT F and T f Ca, respectively, independently of the underlying gauge symmetry. 

Three-jet cross sections were previously measured in ^p [5] and in NC DIS [6,7]. The 
shape of the measured cross sections was well reproduced by perturbative QCD (pQCD) 
calculations and a value of a s was extracted [6]. In this paper, measurements of angu- 
lar correlations in three-jet events in 7£> and NC DIS are presented. The comparison 
between the measurements and fixed-order 0(a 2 ) and 0(a^) perturbative calculations 
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based on different colour configurations provides a stringent test of pQCD predictions 
directly beyond LO and gives insight into the underlying group symmetry. Phase-space 
regions where the angular correlations show potential sensitivity to the presence of the 
TGV were identified. 

2 Theoretical framework 

The dynamics of a gauge theory such as QCD are completely defined by the commutation 
relations between its group generators T\ 

[T\T j ] = iJ2f jk -T k , 

k 

where f l ^ k are the structure constants. The generators T l can be represented as matrices. 
In perturbative calculations, the average (sum) over all possible colour configurations in 
the initial (final) states leads to the appearance of combinatoric factors Cp, Ca and T F , 
which are defined by the relations 

Y^t'A = s ^c F , j2f jkm f jkn = S mn C A , 

k,rj j,k 

Erpm rpn zmnrp 
J-a^Pa — 6 1 F- 

Measurements of the ratios between the colour factors allow the experimental determi- 
nation of the underlying gauge symmetry of the strong interactions. For SU(iV), the 
predicted values of the colour factors are: 

N 2 - 1 

C A = N, C F = and T F = 1/2, 

where N is the number of colour charges. In particular, SU(3) predicts Ca/C f = 9/4 
and T F /C F = 3/8. In contrast, an Abelian gluon theory based on U(l) 3 would predict 
Ca/Cf = and Tf/Cf = 3. A non-Abelian theory based on SO(3) predicts Ca/Cf = 1 
and Tp/Cp = 1- 

The 0(a 2 ) calculations of three-jet cross sections for direct 7p and NC DIS processes can 
be expressed in terms of Ca, Cf and Tp as [8]: 

Cep^Sjcts = C 2 F ■ a a + CfCa ■ o b + CpT F ■ <J C + T f Ca ■ CTDi (1) 
where a a, <Jd are the partonic cross sections for the different contributions (see Fig. 1). 
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3 Definition of the angular correlations 



Angular-correlation observables were devised to distinguish the contributions from the 
different colour configurations. They are defined in terms of the three jets with highest 
transverse energy in an event and the beam direction as: 

• Oh, the angle between the plane determined by the highest-transverse-energy jet and 
the beam and the plane determined by the two jets with lowest transverse energy [9]; 

• ai23, the angle between the two lowest-transverse-energy jets. This variable is based 
on the angle a^ e for e + e~ — > 4 jets [2]; 

• /3ksw, the angle defined via the equation 

cos(/3 K sw) = cos [\ x p 3 ), (p 2 x p B )\ + Z[(ft x p B ), (p 2 x p 3 )))], 

where Pi, i — 1, 3 is the momentum of jet i and ps is a unit vector in the direction of 
the beam; the jets are ordered according to decreasing transverse energy. This variable 
is based on the Korner-Schierholz-Willrodt angle "^ksw f° r e+e ~ ~~ > 4 jets [10]; 

• ?7^*x' the maximum pseudorapidity of the three jets. 

For three-jet events in ep collisions, the variable &n was designed [9] to be sensitive to the 
TGV in quark-induced processes (see Fig. IB). In e + e~ annihilation into four-jet events, 
the distribution of $ksw ^ s sensitive to the differences between qqgg and qqqq final states 
whereas that of a^ e distinguishes between contributions from double-bremsstrahlung 
diagrams and diagrams involving the TGV. 

4 Experimental set-up 

The data samples used in this analysis were collected with the ZEUS detector at HERA 
and correspond to an integrated luminosity of 44.9±0.8 (65.1 ±1.5) pb" 1 for e + p collisions 
taken during 1995-97 (1999-2000) and 16.7 ± 0.3 pb _1 for e~p collisions taken during 
1998-99. During 1995-97 (1998-2000), HERA operated with protons of energy E p = 820 
(920) GeV and positrons or electrons 1 of energy E e = 27.5 GeV, yielding a centre-of-mass 
energy of ^Ts = 300 (318) GeV. 

A detailed description of the ZEUS detector can be found elsewhere [11, 12]. A brief 
outline of the components that are most relevant for this analysis is given below. Charged 
particles were tracked in the central tracking detector (CTD) [13], which operated in a 
magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD consisted 

1 Here and in the following, the term "electron" denotes generically both the electron (e~) and the 
positron (e + ). 
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of 72 cylindrical drift- chamber layers, organised in nine superlayers covering the polar- 
angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks 
was parameterised as <t(pt)/pt = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. The 
tracking system was used to measure the interaction vertex with a typical resolution along 
(transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the energy scale of 
the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [14] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers 
and longitudinally into one electromagnetic section (EMC) and either one (in RCAL) 
or two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the 
calorimeter was called a cell. Under test-beam conditions, the CAL single-particle relative 
energy resolutions were a(E)/E = 0.18 /\^E for electrons and a(E)/E = 0.35/v^ for 
hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'yp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [15], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

5 Data selection and jet search 

A three-level trigger system was used to select events online [12,16]. At the third level, jets 
were reconstructed using the energies and positions of the CAL cells. Events with at least 
one (two) jet(s) with transverse energy in excess of 10 (6) GeV and pseudorapidity below 
2.5 were accepted. For trigger-efficiency studies, no jet algorithm was applied and events 
with a total transverse energy, excluding the energy in the eight CAL towers immediately 
surrounding the forward beampipe, of at least 25 GeV were selected in the 7p sample; for 
the NC DIS sample, events were selected in which the scattered-electron candidate was 
identified using localised energy depositions in the CAL. 

In the offline selection, a reconstructed event vertex consistent with the nominal interac- 
tion position was required and cuts based on tracking information were applied to reduce 
the contamination from beam-induced and cosmic-ray background events. The selection 
criteria of the 7p and NC DIS samples were analogous to previous publications [17, 18]. 

The selected jp sample consisted of events from ep interactions with Q 2 < 1 GeV 2 and 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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a median Q 2 « 10 3 GeV 2 . The event sample was restricted to the kinematic range 
0.2 < y < 0.85, where y is the inelasticity. 

Events from NC DIS interactions were selected from the 1998-2000 data. Two samples 
were studied: Q 2 > 125 GeV 2 and 500 < Q 2 < 5000 GeV 2 . For both samples, |cos7/>| 
was restricted to be below 0.65, where jh, which corresponds to the angle of the scattered 
quark in the quark-parton model, is defined as 

cos = (1 ~y)xE p -yE e 
(1 - y)xE p + yE e 

and x is the Bjorken variable. 

The kx cluster algorithm [19] was used in the longitudinally invariant inclusive mode [20] 
to reconstruct jets in the measured hadronic final state from the energy deposits in the 
CAL cells (calorimetric jets). The axis of the jet was defined according to the Snowmass 
convention [21]. 

For 7p events, the jet search was performed in the r] — <f> plane of the laboratory frame. 
Corrections [17] to the jet transverse energy, Ej, , were applied to the calorimetric jets 
as a function of the jet pseudorapidity, rf 6i , and E^ and averaged over the jet azimuthal 
angle. Events with at least three jets of E^ > 14 GeV and —1 < rf ct < 2.5 were retained. 
Direct 7p events were further selected by requiring a:° bs > 0.8, where x° bs , the fraction of 
the photon momentum participating in the production of the three jets with highest E^, 
is defined as 

if = — ( 4 etl e-^' etl + E^e-^ 2 + 4 et3 e-^ 8t3 ) . 
2yE e V / 

The final 7p data sample contained 1888 events. 

For NC DIS events, the fcr jet algorithm was applied after excluding those cells associated 
with the scattered-electron candidate and the search was conducted in the Breit frame. 
Jet transverse-energy corrections were computed using the method developed in a previous 
analysis [18]. Events were required to have at least three jets satisfying -E^b > & GeV, 
4g' 3 > 5 GeV and —2 < 7^* < 1.5, where £jf B and rf^ are the jet transverse energy 
and pseudorapidity in the Breit frame, respectively. The final NC DIS data sample with 
Q 2 > 125 (500 < Q 2 < 5000) GeV 2 contained 1095 (492) events. 

6 Monte Carlo simulation 

Samples of Monte Carlo (MC) events were generated to determine the response of the 
detector to jets of hadrons and the correction factors necessary to obtain the hadron-level 
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jet cross sections. The hadron level is defined by those hadrons with lifetime r > 10 ps. 
For the NC DIS sample, the MC events were also used to correct the measured cross 
sections for QED radiative effects and the running of a em . 

The generated events were passed through the Geant 3.13-based [22] ZEUS detector- 
and trigger-simulation programs [12]. They were reconstructed and analysed by the same 
program chain as the data. The kx jet algorithm was applied to the MC simulated events 
using the CAL cells in the same way as for the data. The jet algorithm was also applied to 
the final-state particles (hadron level) and the partons available after the parton shower 
(parton level). 

The programs Pythia 6.1 [23] and Herwig 6.1 [24] were used to generate jp events for 
resolved and direct processes. Events were generated using GRV-HO [25] for the photon 
and CTEQ4M [26] for the proton parton distribution functions (PDFs). In both genera- 
tors, the partonic processes are simulated using LO matrix elements, with the inclusion 
of initial- and final-state parton showers. Fragmentation into hadrons is performed using 
the Lund string model [27] as implemented in Jetset [23,28] in the case of Pythia, and 
a cluster model [29] in the case of Herwig. 

Neutral current DIS events including radiative effects were simulated using the Her- 
acles 4.6.1 [30] program with the Djangoh 1.1 [31] interface to the hadronisation 
programs. HERACLES includes corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade is simulated using 
the colour-dipole model (CDM) [32] including the LO QCD diagrams as implemented in 
Ariadne 4.08 [33]; additional samples were generated with the MEPS model of Lepto 
6.5 [34]. Both MC programs use the Lund string model for the hadronisation. The 
CTEQ5D [35] proton PDFs were used for these simulations. 



7 Fixed-order calculations 

The calculations of direct jp processes used in this analysis are based on the program 
by Klasen, Kleinwort and Kramer (KKK) [36]. The number of flavours was set to five; 
the renormalisation, /i R , and factorisation scales, /ip, were set to Hr = = -E™ ax , 
where E™ 3 * is the highest in an event. The calculations were performed using the 
ZEUS-S [37] parameterisations of the proton PDFs; a s was calculated at two loops using 
A^ = 226 MeV, which corresponds to a s (Mz) = 0.118. These calculations are 0{ofy 
and represent the lowest-order contribution to three-jet jp. Full 0(a^) corrections are 
not yet available for three-jet cross sections in jp. 

The calculations of NC DIS processes used in this analysis are based on the program 
NlojetH — h [38], which provides 0(a 2 s ) and 0(a^) predictions for three-jet cross sections. 



6 



The scales were chosen to be /1r = /ip = Q. Other parameters were set as for the 7p 
program. 

In general, the programs mentioned above are very flexible and provide observable- 
independent computations that allow a complete analytical cancellation of the soft and 
collinear singularities encountered in the calculations of jet cross sections. However, these 
programs were written assuming the SU(3) gauge group and the different ingredients 
necessary to perform a calculation according to Eq. (1) were not readily available. The 
programs were rewritten in order to disentagle the colour components to make separate 
predictions for a a, &d- 

The jet algorithm was applied to the partons in the events generated by KKK and 
NlojetH — [- in order to compute the jet cross-section predictions. Thus, these predictions 
refer to jets of partons. Since the measurements refer to jets of hadrons, the calculations 
were corrected to the hadron level. The multiplicative correction factors, defined as the 
ratios between the cross section for jets of hadrons and that for jets of partons, were 
estimated using the MC samples described in Section 6. The normalised cross-section 
calculations changed typically by less than ±5 (10)% for the predictions in jp (NC DIS) 
upon application of the parton-to-hadron corrections. Therefore, the effect of the parton- 
to-hadron corrections on the angular distributions is small. In NC DIS processes, other 
effects not accounted for in the calculations, namely Z° exchange, were also corrected for 
using the MC samples. 

The predictions for jet cross sections are expressed as the convolution of the PDFs and 
the matrix elements, which depend on a s . Both the PDFs and a s evolve with the energy 
scale. In the calculations performed for this analysis, QCD evolution via the DGLAP and 
the renormalisation group equations, respectively, were used. These evolution equations 
also depend on the colour factors. This procedure introduces an additional dependence 
on the colour factors with respect to that shown in Eq. (1); this dependence is suppressed 
by considering normalised cross sections (see Section 8 for the definition of the cross 
sections). The remaining dependence was estimated by comparing to calculations with 
fixed /jj? or /ir. The values chosen for fip and fiR were the mean values of the data 
distributions, (£^ ax ) da ta = 27.8 GeV for 7 p and ^(Q 2 )^ = 31.3 (36.6) GeV for NC DIS 
with Q 2 > 125 (500 < Q 2 < 5000) GeV 2 . 

Figure 2 shows the relative difference of the 0(a 2 ) 'jp calculations with fi F (fi R ) fixed 3 
to those in which /ip = _E™ ax (/i^ = £j? ax ) as a function of the angular variables stud- 
ied. Figures 3(a) to 3(d) show the same relative difference for the 0(a 2 ) NlojetH — h 
calculations for Q 2 > 125 GeV 2 . 

Very small differences are observed for the variation. Sizeable differences for the /i# 
3 When [ip was fixed, [Ir was allowed to vary with the scale, and vice- versa. 
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variation are seen in some regions; in particular, a trend is observed for the relative 
difference as a function of r]^ ax : this trend is due to the fact that the mean values of Q 2 
in each bin of 77^^ increase as r]^ ax decreases. 

These studies demonstrate that the normalised cross sections have little sensitivity to 
the evolution of the PDFs. However, there is still some sensitivity to the running of 
a s . Figures 3(e) to 3(h) show the relative difference for 500 < Q 2 < 5000 GeV 2 . The 
restriction of the phase space further reduces the dependence on the running of a s ; thus, 
this region is more suitable to extract the colour factors in NC DIS at 0(a 2 ). At 0(oP s ) 
(see Fig. 4), the effect due to the running of a s is already very small for Q 2 > 125 GeV 2 . 
Therefore, the wider phase-space region can be kept in an extraction of the colour factors 
at G(a 3 s ). 

The following theoretical uncertainties were considered (as an example of the size of the 
uncertainties, an average value of the effect of each uncertainty on the normalised cross 
section as a function of 9h is shown in parentheses for 7p, NC DIS with Q 2 > 125 GeV 2 
and NC DIS with 500 < Q 2 < 5000 GeV 2 ): 

• the uncertainty in the modelling of the parton shower was estimated by using different 
models (see Section 6) to calculate the parton-to-hadron correction factors (±2.8%, 
±2.9% and ±5.8%); 

• the uncertainty on the calculations due to higher-order terms was estimated by varying 

by a factor of two up and down (11q^%, ±1.6% and ±2.2%); 

• the uncertainty on the calculations due to those on the proton PDFs was estimated 
by repeating the calculations using 22 additional sets from the ZEUS analysis [37]; 
this analysis takes into account the statistical and correlated systematic experimental 
uncertainties of each data set used in the determination of the proton PDFs (±0.7%, 
±0.2% and ±0.1%); 

• the uncertainty on the calculations due to that on a s [Mz) was estimated by repeating 
the calculations using two additional sets of proton PDFs, for which different values 
of a s [Mz) were assumed in the fits. The difference between the calculations using 
these various sets was scaled to reflect the uncertainty on the current world average 
of a s [39] (negligible in all cases); 

• the uncertainty of the calculations due to the choice of /ip was estimated by varying 
/if by a factor of two up and down (negligible in all cases). 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. The dominant source of theoretical uncertainty is that on the 
modelling of the parton shower. 
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8 Definition of the cross sections 



Normalised differential three-jet cross sections were measured as functions of 6h, «23 and 
Pksw using the selected data samples in 'jp and NC DIS. For NC DIS, the normalised 
differential three-jet cross section as a function of ?7^ x was also measured. The normalised 
differential three-jet cross section in bin i for an observable A was obtained using 



1 dOj _ 1 iVdata,i N MC,i 

a dA~ a £■ AAi' N$& fi ' 

where iVdata,; is the number of data events in bin i, (N^^ is the number of MC 

events at hadron (detector) level, C is the integrated luminosity and AAj is the bin width. 
The integrated three-jet cross section, a, was computed using the formula: 



E-^data.i ^MC,i 
7 ' Ardet ' 
j *- JV MC,i 

where the sum runs over all bins. 

For the ^p sample, due to the different centre-of-mass energies of the two data sets used in 
the analysis, the measured normalised differential three-jet cross sections were combined 
using 

^comb _ °"300 " £300 + P"318 " £318 
£300 + £318 

where is the luminosity and a^q is the measured cross section corresponding to 
\/s = 300 or 318 GeV. This formula was applied for combining the differential and total 
cross sections. The same formula was used for computing the 0(a 2 s ) predictions in jp. 



9 Acceptance corrections and experimental uncer- 
tainties 



The Pythia (MEPS) MC samples were used to compute the acceptance corrections to 
the angular distributions of the 7^ (NC DIS) data. These correction factors took into 
account the efficiency of the trigger, the selection criteria and the purity and efficiency 
of the jet reconstruction. The samples of Herwig and CDM were used to compute the 
systematic uncertainties coming from the fragmentation and parton-shower models in 
and NC DIS, respectively. 
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The data E^ , if et and x° bs distributions of the sample, before the x° bs > 0.8 require- 
ment, are shown in Fig. 5 together with the MC simulations of Pythia and Herwig. 
Considering that three-jet events in the MC arise only from the parton-shower approx- 
imation, the description of the data is reasonable. Figure 5(d) shows the resolved and 
direct contributions for the Pythia MC separately. It is observed that the region of 
x° bs > 0.8 is dominated by direct 7p events. The remaining contribution in this region 
from resolved-photon events was estimated using Pythia (Herwig) simulated events to 
be « 25 (31)%. 

Figure 6 shows the data distributions as functions of 9h, «23 and /3ksw together with the 
simulations of Pythia and Herwig for x° bs > 0.8. The Pythia MC predictions describe 
the data distributions well, whereas the description given by Herwig is somewhat poorer. 
It was checked that the angular distributions of the events from resolved processes with 
x obs > o.8 were similar to those from direct processes (see Fig. 7) and, therefore, no 
subtraction of the resolved processes was performed when comparing to the fixed-order 
calculations described in Section 7. 

The data By* B , E^' 3 , rf^ and Q 2 distributions of the NC DIS samples are shown in 
Fig. 8 (9) for Q 2 > 125 (500 < Q 2 < 5000) GeV 2 together with the MC simulations from 
the MEPS and CDM models. Both models give a reasonably good description of the data 
in both kinematic regions. The data distributions of 9h, ^23; Pksw an d rfmax are shown in 
Fig. 10 (11) for Q 2 > 125 (500 < Q 2 < 5000) GeV 2 . The MEPS MC predictions describe 
the data distributions well, whereas the description given by CDM is somewhat poorer. 

A detailed study of the sources contributing to the experimental uncertainties was per- 
formed [40]. The following experimental uncertainties were considered for 7p (as an 
example of the size of the uncertainties, an average value of the effect of each uncertainty 
on the cross section as a function of 9h is shown in parentheses): 

• the effect of the modelling of the parton shower and hadronisation was estimated by 
using Herwig instead of Pythia to evaluate the correction factors (±6.1%); 

• the effect of the uncertainty on the absolute energy scale of the calorimetric jets was 
estimated by varying E^ in simulated events by its uncertainty of ±1%. The method 
used was the same as in earlier publications [17,18,41] (±1.6%); 

• the effect of the uncertainty on the reconstruction of y was estimated by varying its 
value in simulated events by the estimated uncertainty of ±1% (±1.0%); 

• the effect of the uncertainty on the parameterisations of the proton and photon PDFs 
was estimated by using alternative sets of PDFs in the MC simulation to calculate the 
correction factors (±0.4% and ±2.0%, respectively); 

• the uncertainty in the cross sections due to that in the simulation of the trigger 
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(±0.4%). 



For NC DIS events, the following experimental uncertainties were considered (as an ex- 
ample of the size of the uncertainties, an average value of the effect of each uncertainty 
on the cross section as a function of Oh is shown in parentheses for the Q 2 > 125 GeV 2 
and 500 < Q 2 < 5000 GeV 2 kinematic regions): 

• the effect of the modelling of the parton shower was estimated by using CDM instead 
of MEPS to evaluate the correction factors (±5.6% and ±9.1%); 

• the effect of the uncertainty on the absolute energy scale of the calorimetric jets was 
estimated by varying in simulated events by its uncertainty of ±1% for E^ > 
10 GeV and ±3% for lower E^ values (±2.3% and ±1.7%); 

• the uncertainties due to the selection cuts was estimated by varying the values of the 
cuts within the resolution of each variable (less than ±1.6% and less than ±4.2% in 
all cases); 

• the uncertainty on the reconstruction of the boost to the Breit frame was estimated 
by using the direction of the track associated with the scattered electron instead of 
that derived from the impact position as determined from the energy depositions in 
the CAL (±1.6% and ±1.6%); 

• the uncertainty in the absolute energy scale of the electron candidate was estimated 
to be ±1% [42] (±0.2% and ±0.3%); 

• the uncertainty in the cross sections due to that in the simulation of the trigger (±0.5% 
and ±0.5%). 

The effect of these uncertainties on the normalised differential three-jet cross sections is 
small compared to the statistical uncertainties for the measurements presented in Sec- 
tion 10. The systematic uncertainties were added in quadrature to the statistical uncer- 
tainties. 

10 Results 

Normalised differential three-jet cross sections were measured in 7p in the kinematic region 
Q 2 < 1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8. The cross sections were determined for 
jets of hadrons with E^ > 14 GeV and -1 < r/ jct < 2.5. In NC DIS, the cross sections 
were measured in two kinematic regimes: Q 2 > 125 GeV 2 and 500 < Q 2 < 5000 GeV 2 . 
In both cases, it was required that | cos 7^1 < 0.65. The cross sections correspond to jets 
of hadrons with E^ > 8 GeV, E^ 3 > 5 GeV and -2 < iff < 1.5. 
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10.1 Colour components and the triple-gluon vertex 

Normalised differential three-jet cross sections at 0(a 2 ) of the individual colour compo- 
nents from Eq. (1), a a, <?d, were calculated using the programs described in Section 7 
and are shown separately in Fig. 12 for 7p and in Fig. 13 (14) for NC DIS with Q 2 > 125 
(500 < Q 2 < 5000) GeV 2 as functions of the angular variables. In these and subsequent 
figures, the predictions were obtained by integrating over the same bins as for the data. 
The curves shown join the points and are a result of a cubic spline interpolation, except 
in the case of rf^^, for which adjacent points are connected by straight lines. 

The component which contains the contribution from the TGV in quark-induced pro- 
cesses, erg, has a very different shape than the other components for all the angular 
variables considered. The other components have distributions in /Sksw and 9h that 
are similar and are best separated by the distribution of a 2 3 i n IP- in NC DIS with 
500 < Q 2 < 5000 GeV 2 , the different colour components as functions of 6h and Pksw 
also display different shapes. In particular, the <7d component, which also contains a 
TGV, shows a distinct shape for these distributions. This demonstrates that the three-jet 
angular correlations studied show sensitivity to the different colour components. 

In 7 p (NC DIS: Q 2 > 125 GeV 2 , 500 < Q 2 < 5000 GeV 2 ), the SU(3)-based predictions 
for the relative contribution of each colour component are: (A): 0.13 (0.23, 0.30), (B): 
0.10 (0.13, 0.14), (C): 0.45 (0.39, 0.35) and (D): 0.32 (0.25, 0.21). Therefore, the overall 
contribution from the diagrams that involve a TGV, B and D, amounts to 42 (38, 35)% 
in SU(3). 

10.2 Three-jet cross sections in 7p 

The integrated three-jet cross section in 7p in the kinematic range considered was mea- 
sured to be: 

^sjcts = 14.59 ± 0.34 (stat.) t\f x (syst.) pb. 

The predicted 0(a 2 ) integrated cross section, which is the lowest order for this process 
and contains only direct processes, is 8.90 pb- 

The measured normalised differential three-jet cross sections are presented in Fig. 15 and 
Tables 1 to 3 as functions of 9h, cos(a 2 3) and cos(/3ksw)- The measured cross section 
shows a peak at 9h ~ 60°, increases as cos(«23) increases and shows a broad peak in the 
range of cos(/3ksw) between —0.5 to 0.1. 
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10.3 Three-jet cross sections in NC DIS 

The integrated three-jet cross sections in NC DIS for Q 2 > 125 GeV 2 and 500 < Q 2 < 
5000 GeV 2 were measured to be: 

cTep^sjets = 11-48 ± 0.35 (stat.) ± 1.98 (syst.) pb 

and 

a ep ^ 3jets = 5.73 ± 0.26 (stat.) ± 0.60 (syst.) pb. 

The predicted 0(ot\) integrated cross sections are 14.14 ± 3.40 pb and 6.86 ± 1.77 pb for 
the two kinematic regions, respectively. 

The measured normalised differential three-jet cross sections in NC DIS for Q 2 > 125 GeV 2 
and 500 < Q 2 < 5000 GeV 2 are presented in Figs. 16 and 17, respectively, as functions of 
H , cos(a 2 3), cos(/?ksw) an d ( see Tables 4 to 7). The measured cross sections have 
similar shapes in the two kinematic regions considered, except for the distribution as a 
function of cos(/3ksw) : the cross section decreases as cos(/?ksw) increases for 500 < Q 2 < 
5000 GeV 2 whereas for Q 2 > 125 GeV 2 it shows an approximately constant behaviour 
for — 1 < cos(/?ksw) < 0.25. The measured cross section as a function of cos(«23) peaks 
around 0.5 and increases as Oh and rfi^ increase. 

10.4 Comparison to fixed-order calculations 

Calculations at 0(a 2 ) in which each colour contribution in Eq. (1) was weighted accord- 
ing to the colour factors predicted by SU(3) (Cf = 4/3, Ca = 3 and Tp = 1/2) are 
compared to the measurements in Figs. 15 to 19. The theoretical uncertainties are shown 
in Figs. 15, 18 and 19 as hatched bands. Since the calculations are normalised to unity, 
the uncertainties are correlated among the points; this correlation is partially responsible 
for the pulsating pattern exhibited by the theoretical uncertainties. The predictions based 
on SU(3) give a reasonable description of the data for all angular correlations. For 7p, the 
predictions do not include resolved processes (see Section 7), as calculations separated 
according to the different colour factors are not available. Monte Carlo simulations of 
such processes show that their contribution is most likely to be different from that of 
direct processes in the fifth and last bin of (l/a)(da/dcos(a23)) (see Figs. 7b and 15b). 

To illustrate the sensitivity of the measurements to the colour factors, calculations based 
on different symmetry groups are also compared to the data in Figs. 15 to 17. In these 
calculations, the colour components were combined in such a way as to reproduce the 
colour structure of a theory based on the non-Abelian group SU(iV) in the limit of large 
N (C F = 1, C A = 2 and T F = 0), the Abelian group U(l) 3 [C F = 1, C A = and T F = 3), 
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the non-Abelian group SO(3) {Cp = 1/3, Ca = 3 and T F = 1/3) and, as an extreme 
choice, a calculation with Cp = 0. The shapes of the distributions predicted by U(l) 3 
in 7£> are very similar to those by SU(3) due to the smallness of the component Ob and 
the difficulty to distinguish the component op>- In NC DIS, the predictions of U(l) 3 show 
differences of around 10% with respect to those of SU(3), which are of the same order as 
the statistical uncertainties. In both regimes, the data clearly disfavour a theory based 
on SU(iV) in the limit of large N or on Cp = 0. 

Figures 18 and 19 show the measurements in NC DIS compared to the predictions of 
QCD at 0(a 2 ) and 0(a^). This comparison provides a very stringent test of pQCD. The 
0{a1) calculations give a very good description of the data. In particular, a significant 
improvement in the description of the data can be observed for the first bin of the 023 
distribution (Figs. 18b and 19b). 



11 Summary and conclusions 

Measurements of angular correlations in three-jet 7/) and NC DIS were performed in 
ep collisions at HERA using 127 pb -1 of data collected with the ZEUS detector. The 
cross sections refer to jets identified with the kp cluster algorithm in the longitudinally 
invariant inclusive mode and selected with > 14 GeV and — 1 < rf ct < 2.5 (jp) and 
.E^b > 8 GeV, E^f 3 > 5 GeV and -2 < rf^ < 1.5 (NC DIS). The measurements were 
made in the kinematic regions defined by Q 2 < 1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8 (7^0) 
and Q 2 > 125 GeV 2 or 500 < Q 2 < 5000 GeV 2 and | cos7 & | < 0.65 (NC DIS). Normalised 
differential three-jet cross sections were measured as functions of 0h, «23, Pksw and r]^ a ^. 

The colour configuration of the strong interaction was studied for the first time in ep colli- 
sions using the angular correlations in three-jet events. While the extraction of the colour 
factors will require the full analysis of all HERA data and complete 0(al) calculations, 
the studies presented in this paper demonstrate the potential of the method. 

Fixed-order calculations separated according to the colour configurations were used to 
study the sensitivity of the angular correlations to the underlying gauge structure. The 
predicted distributions of 6h, «23 and Pksw clearly isolate the contribution from the 
triple-gluon coupling in quark-induced processes while rf^ x isolates the contribution from 
gluon-induced processes. The variable 023 provides additional separation for the other 
contributions. Furthermore, the studies performed demonstrate that normalised cross 
sections in three-jet ep collisions have reduced sensitivity to the assumed evolution of the 
PDFs and the running of a s . 

The data clearly disfavour theories based on SU(iV) in the limit of large N or Cp = 0. 
Differences between SU(3) and U(l) 3 are smaller than the current statistical uncertainties. 
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The measurements are found to be consistent with the admixture of colour configurations 
as predicted by SU(3). The O(a^) calculations give a very good description of the NC 
DIS data. 
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9 H bin (deg) 


(l/(j) da/d9 H 


"stat 


"syst 


Chad 


0, 9 


0.00264 


0.00038 


±0.00052 


0.93 


9, 18 


0.00393 


0.00044 


±0.00021 


0.94 


18, 27 


0.00507 


0.00051 


+0.00040 
— U.UUUoy 


1.00 


27, 36 


0.00838 


0.00064 


+0.00105 
-0.00104 


0.93 


36, 45 


0.01071 


0.00075 


±0.00023 


0.96 


45, 54 


0.01486 


0.00087 


+0.00021 
-0.00016 


0.94 


54, 63 


0.01795 


0.00098 


+0.00036 
-0.00035 


0.95 


63, 72 


0.01765 


0.00095 


±0.00062 


0.94 


72, 81 


0.01517 


0.00088 


+0.00081 
-0.00084 


0.94 


81, 90 


0.01473 


0.00086 


+0.00075 
-0.00077 


0.96 



Table 1: Normalised differential ep cross section for three-jet photoproduction 
integrated over > 14 GeV and —1 < rf et < 2.5 in the kinematic region defined 
byQ 2 < 1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8 as a function of 9h ■ The statistical 
and systematic uncertainties are shown separately. The multiplicative corrections 
for hadronisation effects to be applied to the parton-level QCD differential cross 
section, dad? o,re shown in the last column. 



cos(tt23) bin 


(1/cr) dcr/(icos(a23) 


<5stat 


^syst 


Chad 


-1, -0.8 


0.0138 


0.0046 


±0.00042 


1.04 


-0.8, -0.6 


0.078 


0.012 


+0.004 
-0.003 


0.96 


-0.6, -0.4 


0.198 


0.022 


+0.026 
-0.027 


0.95 


-0.4, -0.2 


0.343 


0.029 


+0.041 
-0.040 


0.93 


-0.2, 


0.360 


0.029 


±0.010 


0.97 


0, 0.2 


0.512 


0.034 


+0.014 
-0.013 


0.98 


0.2, 0.4 


0.618 


0.037 


+0.015 
-0.016 


1.00 


0.4, 0.6 


0.847 


0.044 


±0.013 


0.99 


0.6, 0.8 


0.937 


0.045 


+0.043 
-0.042 


0.99 


0.8, 1 


1.092 


0.049 


+0.019 
-0.018 


1.02 



Table 2: Normalised differential ep cross section for three-jet photoproduction 
integrated over E ] f l > 14 GeV and — 1 < jf ct < 2.5 in the kinematic region defined 
by Q 2 < 1 GeV 2 , 0.2 < y < 0.85 and x° hs > 0.8 as a function o/cos(a 2 3)- Other 
details as in the caption to Table 1. 
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cos(/3ksw) bin 


(1/cr) da Id cos(/?ksw) 




syst 


Chad 


-1 -0 8 


552 






Q7 


u.o ; u.u 








n qq 


-0 6 -0 4 


745 


042 


+0.032 
-0.031 


97 


-0.4, -0.2 


0.741 


0.042 


±0.039 


0.93 


-0.2, 


0.784 


0.042 


+0.014 
-0.016 


0.96 


0, 0.2 


0.768 


0.042 


±0.046 


0.95 


0.2, 0.4 


0.500 


0.034 


±0.005 


0.94 


0.4, 0.6 


0.200 


0.022 


±0.021 


0.95 


0.6, 0.8 


0.056 


0.010 


+0.010 

-0.009 


0.85 


0.8, 1 


0.0029 


0.0015 


±0.0037 


0.74 



Table 3: Normalised differential ep cross section for three-jet photoproduction 
integrated over > 14 GeV and —1 < rf ct < 2.5 in the kinematic region defined 
by Q 2 < 1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8 as a function of cos(/3 K sw) ■ Other 
details as in the caption to Table 1. 



9 H bin (deg) 


(1/cx) da/d9 H 


"stat 


^syst 


Cqed 


Chad 


Q 2 > 125 GeV 2 


0, 


18 


0.00372 


0.00046 


±0.00031 


0.92 


0.89 


18, 


36 


0.00770 


0.00056 


±0.00095 


0.88 


0.90 


36, 


54 


0.01291 


0.00072 


±0.00045 


0.96 


0.84 


54, 


72 


0.01438 


0.00074 


±0.00042 


1.00 


0.84 


72, 


90 


0.01686 


0.00077 


±0.00160 


0.99 


0.84 






500 < Q 2 


< 5000 GeV 2 






0, 


18 


0.00481 


0.00076 


±0.00048 


0.88 


0.92 


18, 


36 


0.00993 


0.00094 


±0.00231 


0.95 


0.96 


36, 


54 


0.0141 


0.0011 


±0.0004 


0.92 


0.97 


54, 


72 


0.0134 


0.0011 


±0.0008 


1.03 


0.89 


72, 


90 


0.0133 


0.0011 


±0.0023 


0.96 


0.94 



Table 4: Normalised differential ep cross section for three- jet production in NC 
DIS integrated over E^g > 8 GeV, Elf^ > 5 GeV and —2 < 7/3* < 1.5 in the 
kinematic region given by |cos7h| < 0.65 and Q 2 > 125 GeV 2 or 500 < Q 2 < 
5000 GeV 2 as a function of 9h. The multiplicative corrections applied to the dif- 
ferential measured cross section to correct for QED radiative effects, Cqed ; is also 
shown. The multiplicative corrections for hadronisation effects and the Z° '-exchange 
contribution to be applied to the parton-level QCD differential cross section, Chad, 
are shown in the last column. Other details as in the caption to Table 1. 
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cos(a23) bin 


(1/cr) dcr/d cos(a 2 3) <W 


^syst 


Cqed 


Chad 


Q 2 > 125 GeV 2 


-1, -0.6 


0.117 


0.015 


±0.025 


0.96 


0.90 


-0.6, -0.2 


0.338 


0.028 


±0.035 


1.01 


0.70 


-0.2, 0.2 


0.568 


0.032 


±0.018 


0.90 


0.78 


0.2, 0.6 


0.993 


0.037 


±0.021 


0.95 


0.88 


0.6, 1 


0.484 


0.030 


±0.020 


1.02 


1.01 




500 < Q 2 


< 5000 GeV 2 






-1, -0.6 


0.199 


0.030 


±0.018 


1.04 


0.83 


-0.6, -0.2 


0.381 


0.043 


±0.041 


0.97 


0.75 


-0.2, 0.2 


0.589 


0.047 


±0.074 


0.92 


0.83 


0.2, 0.6 


1.018 


0.055 


±0.061 


0.95 


1.07 


0.6, 1 


0.313 


0.036 


±0.022 


0.97 


1.16 



Table 5: Normalised differential ep cross section for three-jet production in NC 
DIS integrated over 

E?tb > 8 GeV, E^' 3 > 5 GeV and -2 < 77-3 < 1.5 in the 
kinematic region given by |cos7/J < 0.65 and Q 2 > 125 GeV 2 or 500 < Q 2 < 
5000 GeV 2 as a function of cos(ct23)- Other details as in the caption to Table 4- 



cos(/3 K sw) bin 


(1/cr) dcr/dcos(j3 


KSw) <Wt 


"syst 


Cqed 


Chad 




Q 2 


> 125 GeV 2 








-1, -0.6 


0.585 




0.031 
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0.92 


0.95 
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0.034 


±0.094 


0.99 
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1.01 


0.85 


0.2, 0.6 


0.332 




0.026 


±0.025 


0.92 


0.74 


0.6, 1 
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0.020 


±0.022 


0.93 


0.71 




500 < Q 2 


< 5000 GeV 2 








-1, -0.6 


0.770 




0.052 


±0.076 


0.94 


1.04 


-0.6, -0.2 


0.536 




0.045 


±0.112 


0.93 


0.97 


-0.2, 0.2 


0.497 




0.045 


±0.037 


1.01 


0.94 


0.2, 0.6 


0.430 




0.044 


±0.058 


1.01 


0.84 


0.6, 1 


0.267 




0.036 


±0.061 


0.89 


0.78 



Table 6: Normalised differential ep cross section for three- jet production in NC 
DIS integrated over -E^b > 8 GeV, E 3 ^' 3 > 5 GeV and —2 < rf^ < 1.5 in the 
kinematic region given by |cos7h| < 0.65 and Q 2 > 125 GeV 2 or 500 < Q 2 < 
5000 GeV 2 as a function o/cos(/9ksw)- Other details as in the caption to Table 4- 
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r? J bin 

/max 


(1/a) da/drf^ 


^stat 


^syst 


Cqed 


Chad 


Q 2 > 125 GeV 2 


-2, -0.1 


0.0042 


0.0013 


±0.0006 


1.07 


0.61 


-0.1, 0.3 


0.092 


0.016 


±0.012 


1.17 


0.77 


0.3, 0.7 


0.267 


0.024 


±0.054 


0.96 


0.81 


0.7, 1.1 


0.751 


0.034 


±0.016 


0.93 


0.83 


1.1, 1.5 


1.370 


0.038 


±0.048 


0.96 


0.88 




500 < Q 2 


< 5000 GeV 2 






-2, -0.1 


0.0059 


0.0021 


±0.0022 


1.14 


0.62 


-0.1, 0.3 


0.110 


0.022 


±0.011 


0.96 


0.77 


0.3, 0.7 


0.378 


0.040 


±0.084 


0.96 


0.86 


0.7, 1.1 


0.918 


0.054 


±0.052 


0.93 


0.93 


1.1, 1.5 


1.066 


0.056 


±0.035 


0.98 


1.00 



Table 7: Normalised differential ep cross section for three- jet production in NC 
DIS integrated over -Ej"b > 8 GeV, E^' 3 > 5 GeV and —2 < rf^ < 1.5 in the 
kinematic region given by |cos7h| < 0.65 and Q 2 > 125 GeV 2 or 500 < Q 2 < 
5000 GeV 2 as a function of rf^ x . Other details as in the caption to Table 4- 
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Figure 1: Examples of diagrams for the photoproduction of three-jet events through 
direct-photon processes and in NC DIS three-jet events in each colour configuration: 
(A) double-gluon bremsstrahlung from a quark line; (B) the splitting of a virtual 
gluon into a pair of final- state gluons; (C) the production of a qq pair through the 
exchange of a virtual gluon emitted by an incoming quark; (D) the production of 
a qq pair through the exchange of a virtual gluon arising from the splitting of an 
incoming gluon. 
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Figure 2: Relative difference between the 0(a 2 s ) calculations with lir = 27.8 GeV 
and the calculations with lir = i?™ ax (dots) and between the 0(ct 2 s ) calculations 
with [If = 27.8 GeV and the calculations with \ip — E^ ax (open circles) in jp 
as functions of (a) 9 H , (b) cos(a 2 3) and (c) cos(/3ksw)- These calculations do not 
include corrections for hadronisation effects. 
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Figure 3: Relative difference between the 0{a 2 s ) calculations with fixed lir and 
the calculations with llr = Q (dots) and between the 0(a 2 s ) calculations with fixed 
Li f an d the calculations with /xp = Q (open circles) in NC DIS as functions of 
(a,e) Oh, (b,f) cos(a23), (c,g) cos(/5ksw) an d (d,h) r]^ x . These calculations do not 
include corrections for hadronisation effects. 
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Figure 4: Relative difference between the 0{oP s ) calculations with fixed fin and 
the calculations with /ir = Q (dots) and between the 0(a^) calculations with fixed 
Hp and the calculations with fi F = Q (open circles) in NC DIS as functions of 
(a,e) Oh, (b,f) cos(«23), (c,g) cos(/5ksw) and (d,h) r]^ x . These calculations do not 
include corrections for hadronisation effects. 
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ZEUS 




Figure 5: Detector-level data distributions for three-jet photoproduction (dots) 
with > 14 GeV and -1 < rf ct < 2.5 in the kinematic region given by Q 2 < 
1 GeV 2 and 0.2 < y < 0.85 as functions of (a) , (b) rf ct and (c,d) xf s . 
For comparison, the distributions of the PYTHIA (solid histograms) and HERWIG 
( dashed histograms ) MC models for resolved plus direct processes normalised to the 
data are included. In (d), the contributions for resolved (dotted histogram) and 
direct (dot-dashed histogram) processes from Pythia MC are shown separately. 
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Figure 6: Detector-level data distributions for three-jet photoproduction (dots) 
with Ej, > 14 GeV and — 1 < rf Qt < 2.5 in the kinematic region given by Q 2 < 
1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8 as functions of (a) 8 H , (b) cos(a 23 ) and 
(c) cos(/5ksw)- Other details as in the caption to Fig. 5. 
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Figure 7: Detector-level data distributions for three-jet photoproduction (dots) 
with E}£ > 14 GeV and —1 < rf ct < 2.5 m the kinematic region given by Q 2 < 
1 GeV 2 , 0.2 < y < 0.85 and x° bs > 0.8 as functions of (a) 8 H , (b) cos(a 2 3) and (c) 
cos(/?ksw)- The predictions for resolved (dotted histogram) and direct (dot-dashed 
histogram) processes from the PYTHIA MC normalised separately to the data are 
also shown. Other details as in the caption to Fig. 5. 
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Figure 8: Detector-level data distributions for three-jet production in NC DIS 
(dots) with E^l > 8 GeV, E^> 3 > 5 GeV and -2 < ?7g < 1.5 in the kinematic 

region given by Q 2 > 125 GeV 2 and | cos7h| < 0.65 as functions of (a) ^ et B , (^j 

-^tb' 3 ' ( c ) ^b* anc ^ W -^ 9r comparison, the distributions of the MEPS (solid 
histograms) and CDM (dashed histograms) MC models normalised to the data are 
included. 
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Figure 9: Detector-level data distributions for three- jet production in NC DIS 
(dots) with -E^b > 8 GeV } E^™' 3 > 5 GeV and —2 < rf^ 1 < 1.5 m £/ie kinematic 
region given by 500 < Q 2 < 5000 GeF 2 and ICOS7/J < 0.65 as functions of (a) 
Etb> 3 ) -S'tb' 3 an d ( c ) Vrt ■ Other details as in the caption to Fig. 8. 
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Figure 10: Detector-level data distributions for three-jet production in NC DIS 
(dots) with -E^g > 8 GeV, E^' 3 > 5 GeV and —2 < rf-^ < 1.5 m i/ie kinematic 
region given by Q 2 > 125 GeV 2 and |cos7/j| < 0.65 as functions of (a) 9 H , (b) 
cos(«23), (c) cos(/?ksw) and (d) if^^. For comparison, the distributions of the 
MEPS (solid histograms) and CDM (dashed histograms) MC models normalised to 
the data are included. 
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Figure 11: Detector-level data distributions for three-jet production in NC DIS 
(dots) with > 8 GeV, E^' 3 > 5 Ge^ and -2 < 773* < 1.5 in the kinematic 
region given by 500 < Q 2 < 5000 GeV 2 and Icc^^l < 0.65 as functions of (a) 
6 H , (b) cos(a 2 3); (c) cos(/?ksw) an d (d) n^* x . Other details as in the caption to 
Fig. 10. 
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Figure 12: Predicted normalised differential ep cross sections for three-jet direct- 
photon processes at 0(a> 2 ) integrated over > 14 GeV and -1 < rf ct < 2.5 in 
the kinematic region defined by Q 2 < 1 GeV 2 and 0.2 < y < 0.85 as functions 
of (a) 9 H , (b) cos(a 2 3) and (c) cos(/9ksw)- I n eac h figure, the predictions for the 
colour components are shown: a a (dashed lines), Ob (solid lines), Oc (dot-dashed 
lines) and o~d (dotted lines). These calculations do not include corrections for 
hadronisation effects. 
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Figure 13: Predicted normalised differential ep cross sections for three-jet pro- 
duction in NC DIS at 0{a 2 s ) integrated over > 8 GeV, E$ t2 ' 3 > 5 GeV 

and —2 < rj^ < 1.5 in the kinematic region given by Q 2 > 125 GeV 2 and 
| cos7h| < 0.65 as functions of (a) 9 H , (b) cos(a 23 ), (c) cos(/5ksw) an d (d) r]^ ax . 
Other details as in the caption to Fig. 12. These calculations do not include cor- 
rections for hadronisation effects. 
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Figure 14: Predicted normalised differential ep cross sections for three-jet pro- 
duction in NC DIS at 0(a 2 ) integrated over £y g > 8 GeV, E$ t2 > 3 > 5 GeV 

and —2 < rf^ < 1.5 in the kinematic region given by 500 < Q 2 < 5000 GeV 2 
and |cos7/J < 0.65 as functions of (a) 9 H , (b) cos(a 2 3) ; (c) cos(/?ksw) an d (d) 
rfmax- Other details as in the caption to Fig. 12. These calculations do not include 
corrections for hadronisation effects. 
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Figure 15: Measured normalised differential ep cross sections for three-jet pho- 
toproduction (dots) integrated over > 14 GeV and -1 < ^ ct < 2.5 in the 
kinematic region defined by Q 2 < 1 GeV 2 , 0.2 < y < 0.85 and a;° bs > 0.8 as func- 
tions of (a) 9 H , (b) cos(a23) and (c) cos(/3ksw)- The data points are plotted at the 
bin centres. The inner error bars represent the statistical uncertainties of the data, 
and the outer error bars show the statistical and systematic uncertainties added 
in quadrature. For comparison, the 0(a> 2 ) calculations for direct-photon processes 
based on SU(3) (solid lines), U(l) 3 (dashed lines), SU(N) in the limit of large 
N (dot-dashed lines), Cp = (short-spaced dotted lines) and SO(3) (long-spaced 
dotted lines) are included. The lower part of the figures displays the relative dif- 
ference to the calculations based on SU(3) and the hatched band shows the relative 
uncertainty of this calculation. 
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Figure 16: Measured normalised differential ep cross sections for three-jet 
production in NC DIS (dots) integrated over E^ > 8 GeV, E$ t2 ' 3 > 5 GeV 

and —2 < rf^ < 1.5 in the kinematic region given by Q 2 > 125 GeV 2 and 
|cos7h| < 0.65 as functions of (a) 6h, (b) cos(a23) ; (c) cos(/?ksw) and (d) rft^. 
Other details as in the caption to Fig. 15. 



38 



ZEUS 



SU(3) 

\3 



• ZEUS DIS 82 pb" 1 
U(ir (a) 
SU(N), large N 

c F = o 



SO(3) 




20 



40 



60 80 

H (deg) 



V 

© 
U 

73 



0.5 



GO 

o 
-3 



o 

0.5 


-0.5 
-1 



500 < Q 2 < 5000 GeV 2 



(b) 




-0.5 



0.5 1 

cos(cc 23 ) 



(c) 




4- es 

-a 



-0.5 



0.5 1 

cos((3 KSW ) 



£> 
o 

"3 
- 



1.5 



0.5 


0.5 


-0.5 



(d) 




: — 1 — 1 — 1 — 1 — * 


■fT" — i — i — — i — i — i — i — — i — h 


- 





max 



Figure 17: Measured normalised differential ep cross sections for three-jet pro- 
duction in NC DIS (dots) integrated over Eifv > 8 GeV, E^ t2 ' 3 > 5 GeV and 

< 1.5 in the kinematic region given by 500 < Q 2 < 5000 GeV 2 and 



-2 < 



ct 



IC0S7/J < 0.65 as functions of (a) 6h, (b) cos(a23) 7 (c) cos(/3ksw) 
Other details as in the caption to Fig. 15. 
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Figure 18: Measured normalised differential ep cross sections for three-jet 
production in NC DIS (dots) integrated over > 8 GeV, E^ t2 ' 3 > 5 GeV 

and —2 < rf^ < 1.5 in the kinematic region given by Q 2 > 125 GeV 2 and 
| cos7h| < 0.65 as functions of (a) 9h, (b) cos^s)? (c) cos(/5ksw) and (d) r/^ ax . 
For comparison, the 0(a 2 ) (dashed lines) and 0(a 3 ) (solid lines) QCD calculations 
are also included. The hatched band displays the relative theoretical uncertainty of 
the 0(a 3 ) calculation. Other details as in the caption to Fig. 15. 



40 



ZEUS 



O(o# 

0(0$ 



• ZEUS DIS 82 pb" 1 



(a) 




I7//7/I theoretical uncertainty 

i I gggj 



20 



40 



60 80 

H (deg) 



o 
u 

e 
-a 



0.5 



^ 

o °- 5 

o 





-0.5 
-1 



500 < Q 2 < 5000 GeV 2 



(b) 




-0.5 



0.5 1 

cos(cc 23 ) 



(c) 





M 



-a 



-0.5 



T3 



0.5 1 

cos((3 KSW ) 



1.5 



0.5 







V 

o 05 

o 

* o 

'■3 -0.5 



^^^^^ 


(d): 


1 1 1 1 


It I I 1 I I 1 1 1 


1 







1 ^jet 
max 



Figure 19: Measured normalised differential ep cross sections for three-jet pro- 
duction in NC DIS (dots) integrated over > 8 GeV, E^ t2 ' 3 > 5 GeV and 

— 2 < T]^* < 1.5 in the kinematic region given by 500 < Q 2 < 5000 GeV 2 and 
| cos7h| < 0.65 as functions of (a) Oh, (b) cos(a;23), (c) cos(/5ksw) and (d) r]^ ax . 
For comparison, the 0(a 2 ) (dashed lines) and 0(a^) (solid lines) QCD calculations 
are also included. The hatched band displays the relative theoretical uncertainty of 
the 0{a z s ) calculation. Other details as in the caption to Fig. 15. 
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